function g=find_g_kumagai(x)
h=-0.365;
c=[0.20173476 730418.72 1000000.0 1.0 26.0];
c1=c(1);
c2=c(2);
c3=c(3);
c4=c(4);
c5=c(5);
x2=(x-h).^2;
g=c1+c2*x2./(c3+x2).*(1+c4*exp(-c5*x2));
